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O . ABSTRACT 

(N ■ 

D ' We develop a model for the infrared, optical and soft X-ray emission of the Crab 

'' . pulsar in terms of anisotropic synchrotron emission by relativistic particles in an outer 

^^ . gap scenario with a single energy distribution A^(7) oc 7~^. It is shown that such a 

\ distribution is naturally produced in an efficient pair cascade and that the energy of the 

'T^ \ primary particles are limited by synchrotron radiation to 7 ~ 10^. It is further shown 

'i^ I that this synchrotron model is able to reproduce the spectral shape between the infrared 

^j| ' and soft X-rays and also the corresponding luminosities. In particular, the long standing 

Q^ ■ problem of the rapid spectral decline towards infrared frequencies is understandable as 

^? ■ emission at very small pitch angles from low energy particles with 7 ~ 10^. 
<^ . Finally, we show that the scaling of our synchrotron model explains the observed 

r"| . correlation between the X-ray luminosity and the spin-down luminosity of the neutron 

0^\ star Lx ~ IQ-^Lsd found by Becker & Triimper (1997). 

O' 

'^ I Subject headings: Stars: Pulsars: Individual (PSR 0531-1-21) — Radiation Mechanisms: 

^ ' Nonthermal (Synchrotron Radiation) — Stars: X-Rays 

>'■ 

^: 

» . 1. Introduction 

The Crab pulsar has a continuous spectrum from the optical to X-rays and 7-rays with different 
power-laws (Lyne &: Graham-Smith 1990). The spectral index en, defined as Ii, oc i^~", varies from 
a = 1.1 in the 7-region via 0.7 in the hard X-ray region to 0.5 at soft X-rays (Toor & Seward 
1977) and zero at optical frequencies. In the far infrared region the spectrum is inverted a = —2 
and cuts off sharply towards lower frequencies. This behaviour is not accompanied by dramatic 
pulse profile changes, as one would expect from saturation or self-absorption effects. Self-absorption 
should first influence the peak intensity. Although Penny (1982) claims to have measured about 
20% peak-depression in the near infrared, the situation is ambiguous. Absorption outside the pulsar 
environment would not influence the pulse form, but the nature of such an absorber is not clear. The 
conventionally used optically thin synchrotron radiation theory yields rising spectra if the particle 
spectrum has a low-energy cutoff, with a spectral index of —1/3, if the angle integrated single 
particle emissivity is used and —2/3, if the angular dependence is retained, e.g., when the emission 



is spread over angles becoming comparable to the pulsar emission cone angle. The infrared spectrum 
measured by Middleditch, Pennypacker & Burns (1983) therefore is still awaiting an explanation. 
It is one of the purposes of this paper to show that the theory of optically thin synchrotron radiation 
at very small pitch angles gives a possible solution to this problem. Epstein &: Petrosian (1973) 
applied small pitch angle synchrotron radiation to the Crab optical pulsations, but the detailed 
spectral features were only calculated for AGN. 



2. Outer Gap Scenario and Particle Energy Spectrum 

We suggest that the infrared, optical and ultraviolett emission (and also the X-ray and 7- 
ray emission) can be explained in terms of an outer magnetospheric gap model (Cheng, Ho & 
Rudermann 1986). Models for the Gamma ray emission from outer gaps have been published by 
Romani and Yadigaroglu (1995) and Chiang and Romani (1994). This idea is also supported by 
the different widths of pulse profiles in the radio (inner magnetosphere) and the higher frequency 
range. The magnetic field strength close to the velocity of light cylinder of Crab pulsar is S ~ lO^G. 

In such a scenario particles are accelerated to high Lorentz factors 7^ ~ 10^ by large potential 
drops along the magnetic field lines, usually thought to be limited by curvature radiation. The 
electric field strength in the gap of a pulsar rotating with angular velocity 0, = 2it/P has been 
calculated by Cheng, Ho & Rudermann (1986) to be 

^» = ij:7=«-^'<'»(io^) (iP^j (33^) (1^0)="™"/""- « 

where Re is the curvature radius of the field lines in the gap region and a is the dimension of the 
gap perpendicular to the magnetic field B. The rate at which particles gain energy in such an 
electric field is 

-Pgain = ceEo , (2) 

and the energy loss rate due to curvature radiation is 

Per = 1^1' = 4.6 X lO-^R-'^' ergs-i . (3) 

We now consider a competing process, which can limit the acceleration of particles to high 
energies, namely synchrotron radiation. The energy loss rate of a single particle is then given as 

Psyn = 1.6 X 10-^^7^52 sin^ ^ erg s"^ , (4) 

which is larger than the curvature radiation losses if the pitch angle is 



where sin "J/ ps vp has been used. 

Since we are interested in the acceleration process close to the light cylinder radius Ric, we 
now calculate the accelerating electric field there. Let r and 9 be the radial distance from the 
neutron star center and the angle between the magnetic moment and the radius vector to a point 
in the magnetosphere, respectively. The angle between the magnetic moment and the rotation axis 
is called a. Let further be re = r/R^s, where R^s ~ 10® cm is the radius of the neutron star. The 
radius of the light cylinder Ric is given by QRi^ = c, so that 

In a purely dipolar magnetic field the relation 

sin^ 6 sin^ ^n 

-r = -ir- ^^^ 

holds for every field line, where ^o defines its footpoint at the neutron star surface. The curvature 
radius of a magnetic field line is known to be 

R, = ^^1/2 . (8) 

In an aligned rotator the last closed field line would be characterized hy 9 = 11/2 at the light 
cylinder radius, whereas in the non-aligned case it is given by ^ = tt/2 — a. It is then easily verified 
with eq. (7), that 

sin ^0 = a^ic cos a . (9) 

We therefore find for the curvature radius at the light cylinder 

i?e = ^^^^^ = -^ . (10) 

cos a cos a 
The accelerating electric field near the light cylinder is then 

Eq = — 2" COS a 5 = 5 cos a 5, (11) 

which defines 5 as the ratio of the gap thickness to the light cylinder radius, in the following assumed 
to be of order 0.1. 

One can now estimate the equilibrium Lorentz factor for which energy gains and synchrotron 
losses are equal, as 

These particles are called primary in the following. By inverse Compton scattering they can produce 
7-rays with energies up to 5 TeV. We use a pitch angle of order ^ ~ 10~^ which in the following 



sections is shown to be required to reproduce the infrared and optical spectrum, as well as the 
luminosity in this band. 

Due to interactions of these 7-rays with photons a pair creation cascade is ignited, which can 
produce several generations of positrons and electrons of equal amount and a certain multiplicity M 
compared to the Goldreich- Julian density. The initial momentum of the created pairs will be parallel 
to that of the primary particle, i.e. the pairs will also have a pitch angle ^ ~ 10~^. In typical cascade 
models the pair density exceeds the density of primaries by a factor of 10^ . . . 10^ (Melrose 1995). 
Numerical simulations of such a cascade showed that a power-law energy distribution N{'y) = Nq^~^ 
will be produced, with a typical index of s ~ 2. The first cascade considerations were performed 
by Daugherty & Harding (1982) for an inner gap scenario but if the cascade is possible in the 
outer gap (which is commonly assumed) the energy spectrum of the secondaries should not be 
too different. The spectrum cuts off below 7 < 10^ due to the energy theshold of the interacting 
photons to produce a pair. 

This particle energy spectrum can also be found by the following simple argument. Consider 
a primary particle with energy ^piTieC^ cascading into M('~f) particles, each with an energy jrUeC^. 
If the cascade is very efficient, energy conservation leads to 

M(7) = ^r^— . (13) 

7 7 

Since M{'y) is the number of particles around Lorentz factor 7, the differential number density, or 
energy spectrum is then 

7V(7) oc ^^ oc ^ . (14) 

7 T 

The initial momentum of the primary particle (parallel to the local magnetic field) defines also 
the momentum of the 7-rays emitted by inverse Compton scattering with ambient photons. The 
7-rays eventually leave the gap because they are not influenced by the magnetic field and proceed 
to created pairs by photon-photon interactions. The subsequent generations of pairs are then no 
longer affected by the electric field of the gap. So it seems reasonable to assume that via creation 
all particles have the same pitch angle with a following evolution of the pitch angle distribution due 
to motion along the magnetic field lines and radiative losses. Of course the emisssion of photons 
will influence the particle energy spectrum, giving a steepening after a loss time scale tsyn, affecting 
first the high energy part. 



3. Synchrotron Emission at Very Small Pitch Angles 

Since the relativistic particles of the pair cascade are created almost in the direction of the 
magnetic field lines one has to be careful to use the correct validity range of the approximations 
made in synchrotron radiation theory. In Epstein (1973) the synchrotron emission process by 



particles with very small pitch angles is discussed in detail. The emissivity at very low pitch angles 
(^ <^ 1/7) of a single particle is given by 
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where ub = eB/2'Kmc = 2.80 x lO^^SeHz is the nonrelativistic gyro frequency and 9 is the angle 
between the magnetic field B {Bq: B in units of 10® G) and the direction of emission. Since we 
restrict our discussion to small pitch and emission angles we set sin^ ~ ^ and sin 9 ~ 6. This 
formula has to be applied when angles 9 < I/7 are resolved in the observations, i.e. in the case of a 
pulsar, when the pulse width A(p becomes comparable to the emission cone angle of the particles, 
1/7 > 27rA0. Since the maximum of the emission is in the forward direction (9 = 0) the emission 
is dominated by the field lines that point towards the observer at each phase of the pulse. The 
degree of circular and linear polarization is given by 
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where ^u = ±1 is the sign of the charge. 
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When the typical angle of emission 9 is small compared to the pulse width it is useful to 
integrate the emissivity over all angles to get the total emission spectrum of a single particle 
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up to z^ < 2yvB- The energy loss rate of a single particle is obtained by integrating the total 
emission over all frequencies 



Lr 



dE 
~dt 



Ey dv 



87r2e2472^2 



'0 3c 

which is also the well known result at large pitch angles, used in eq. (4). 



(19) 



We now calculate the emission from a power law distribution of relativistic particles, iV(7) 
A'^07^'*, in the direction along the magnetic field = according to 
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This result shows that the spectrum rises very steeply, in fact I^ oc v"^ for s = 2. Therefore the 
sharp rise in the spectrum of the Crab pulsar at infrared frequencies can be explained as synchrotron 



emission at very small pitch angles. The Lorentz factors in the case of the Crab then have to be 
of the order of 7 < 1/A(^ rsj lo^ so that this description applies. This is consistent with the 
requirement that 2^ub lies in the infrared part of the spectrum. 

In case of somewhat higher Lorentz factors the single particle beam angle becomes smaller than 
the phase angle of the pulse and the total (or angle averaged) power spectrum has to be applied. 
To find the spectrum from a power-law distribution the monochromatic approximation 

e,^-—6{u-2jUB) (22) 

is used for simplicity. This is a good approximation since the emission is sharply peaked at 2^i'b 
and one may just put all the emission at this frequency. As a result we then find 

U<.J..N(y),-,= Sf^(^) • (23) 

In the case of s = 2 this gives the observed flat, a ~ 0, optical/UV spectrum of the Crab pulsar, 
with 7 rs^ 10^ — 10^. Since the flat spectrum continues up to the far-ultraviolett (Lundqvist et al. 
1999) the pitch angle involved has to be ^ ~ 10~^ so that the emission is still at very low pitch 
angles (7^ < 1). 



4. Crab Pulsar Spectrum from the Infrared to X-rays 

We consider now the injection of a power law distribution of a relativistic electron-positron 
plasma A^(7) = Nq^^"^, which is produced in the cascading process, along the magnetic field lines. 
We interpret the infrared and optical emission of the crab pulsar with its fast rise towards the 
optical as synchrotron emission from relativistic particles with very small pitch angles (7 ^ 1/^)- 
As shown in the previous section we expect a spectral index of a = —2 in the infrared part and 
a = in the optical part, if the pitch angle is of order ^ ~ 10"'^ for Lorentz factors ~ 10^ — 10^. 

Although the circular polarization will be at its maximum when looking along the field lines, 
in the cascade process an equal amount of electrons and positrons will be produced, giving a net 
zero circular polarization. On the other hand the linear polarization becomes maximal for 6 = I/7 
for both species and is zero for 6 = 0. If the emission cone of the pulsar is not cut right through 
the center by the line of sight but slightly offset, then one expects a polarization angle swing. 

The soft X-ray emission in this model comes from the particles with larger pitch angles (^ ^ 
1/7), radiating at a typical frequency 

i^c = li^B^l^ (24) 

The spectrum in this large pitch angle limit is given by 

I, oc ^.-(^-i)/2 (25) 



which yields the observed I^ oc u^^'^ at soft X-rays for s = 2. Therefore a single power-law 
distribution of relativistic electrons and positrons can produce the complete spectrum from infrared 
to X-ray frequencies. 



5. Simple Luminosity Estimates for the Optical and X-ray Emission 

We now estimate the energy radiated by the pulsar in the optical and in the X-ray band simply 
as the number of particles times the energy radiated by one particle. Thus we have 

L = MncjVPsyn (26) 

with 

Psyn = 1.6 X lO-^^^^^^^^ ergs"^ (27) 

being the energy radiated by a single particle. M is the multiplicity of the relativistic particle 
density wihin the volume V as compared to the Goldreich-Julian density 

ncj = 6.9 X IQ-'^BP-^ cm"^ , (28) 

where P is the period of the pulsar in seconds, x = R/R^s is the distance in units of the neutron 
star radius. The magnetic dipole field falls off as i? = Bq/x^, with Bq being the field strength at 
the surface of the neutron star. Since we assume the radiation to occur close to the light cylinder 
we set 

X = xic = 4774 P . (29) 

The volume is estimated as a spherical shell within a fraction / of the light cylinder radius 

V = AttRI ■ fRic = 1.4 X 10^°/^^ cm^ • (30) 

With M = 7p/7 we then find that 

L = 1.5 X 10^^7p7/^2^3p2 g^gg-i _ (3^) 

The frequency of the optical emission (at small pitch angles) is given by 

uopt = 2.8 X lOSopt^ Hz , (32) 

so that, by eliminating 7, we get 

Lop, = 4 X 10- cosV^ a ( ±) f ^ V Al ^ (^) ^ (^^) 

„ X 3/2 / p N -5/2 

' ' ' ergs . (33) 



10i2-5Gy V33ms 

The parameter values chosen for the Crab pulsar are consistent with the considerations of the 
spectral shape and the cascading process, i.e. the particles with 7 ~ 10^ — 10^ have pitch angles 



^opt ~ 10 ^. This crude estimate gives the right order of magnitude of the Crab pulsars optical 
luminosity. 

In the X-ray regime where the large pitch angle approximation for the frequency applies 

i/x = 2.8 X lO^^xTx^ Hz (34) 



we then find 
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ergs . (35) 



In order to get the right X-ray luminosity a pitch angle of order ^x ~ 0.1 is required. To emit at 
soft X-ray frequencies again the Lorentz factors have to be in the range 7 ~ 10^ — 10^. 

This result can now be compared with the spin-down luminosity of the Crab pulsar 

i.<.=3^2xl0"(^)'(^) 'erg,-'. (36) 

The luminosity radiated in X-rays by the synchrotron mechanism shows the same dependence on 
the magnetic field strength and the period as the rotational energy loss of the pulsar. We can 
therefore write 

....x.-»™V.„(_L)(i)(A,)-'(^)-(_J^)-,.,. ,3, 

In order to get the right X-ray luminosity a pitch angle of order ^x ~ 0.1 is required. To emit at 
soft X-ray frequencies again the Lorentz factors have to be in the range 7 ~ 10^ — 10^. 

Since this relation is now independent of Bq and P it should apply to any pulsar, not only 
the Crab pulsar. The observations of X-ray selected pulsars by Becker & Triimper (1997) have 
revealed such a dependence in the form L^ ~ 10~^Lsd- The theoretically found proportionality 
between the X-ray luminosity and the spin-down luminosity strongly supports our assumption, 
that synchrotron emission is the dominant radiation process. On the basis of the model presented 
here, this observation implies that the parameter combination cos^'^ aJ/^p^^x has about the 
same value for every pulsar in their list. 

The luminosity estimates also contain the correct spectral dependencies, because I^ oc L^pi/i' oc 
const, for the optical and I,^ oc L^/v oc u^^''^ for the X-ray emission. 



6. Discussion and Conclusions 

The high frequency emission of the Crab pulsar is a challenge for the theoretical description. 
Earlier investigations put the location of the infrared and optical pulses close to the light cylinder 



radius (Smith et al. 1988) and led us to the idea that these distant regions of the magnetosphere 
may be responsible for the emission up to the X-ray and gamma ray range. First we were able 
to show that within an outer gap description and a standard cascade model we expect a particle 
energy distribution A'^(7) oc 7~^, being the basis for the investigations of the radiative properties 
of the Crab pulsar. The energy of the primary particles accelerated in the electric field of the gap 
is shown to be limited by synchrotron radiation and not curvature radiation to '^ lO^rricC^- This 
holds for pitch angles ^ > 10~^. 

A detailed inspection of an anisotropic synchrotron model led to the result that the emission of 
particles with small pitch angles with this single energy distribution reproduces the key properties 
of the spectrum and luminosity from the infrared up to soft X-rays. This model explains all relevant 
features, including the long-standing problem of the increase of the infrared spectrum I^ oc i^^, the 
flat spectrum I^ oc v^ in the optical and UV range and the X-ray spectrum Ii, oc i'~^'^. The only 
difference concerns the pitch angle of the emitting particles and the (energy dependent) multiplicity 
M = M(7). The infrared/optical/UV part of the spectrum is emitted by particles with Lorentz 
factors 7 -^ 10^ — 10^ and pitch angle ^ ~ 10^^, corresponding to the pitch angle of the primary 
particles accelerated in the electric field of the gap. The X-rays are radiated by particles in the 
same energy range but with a pitch angle \I' ~ 0.1. It is not clear what the fundamental physical 
scenario is, which gives the required values for the pitch angle. 

The X-ray emission from particles with the smaller pitch angles, involving higher energy par- 
ticles, which are smaller in number, would not be seen due to the smaller power radiated. It is 
dominated by particles with the larger pitch angle. On the other hand for pulsars with smaller 
magnetic field strength at the light cylinder, the emission at larger pitch angles might be seen in 
the optical band. In that case the z^~^'^-spectrum is also expected there, then accompanied by a 
luminosity relation similar to that for the X-rays. 

From simple luminosity arguments we find that this model is able to yield the desired energy 
output in the optical as well as in the X-ray part of the spectrum. We also find that the observed 
proportionality of the X-ray flux and the spin-down luminosity is a consequence of the synchrotron 
radiation model for any pulsar that developes an efficient pair cascade. 
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